Disentanglement and decoherence from classical non-Markovian noise: 

Random telegraph noise 
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We calculate the two-qubit disentanglement due to classical random telegraph noise using the 
quasi-Hamiltonian method. This allows us to obtain analytical results even for strong coupling 
and mixed noise, important when the qubits have tunable working point. We determine when 
entanglement sudden death and revival occur as functions of qubit working point, noise coupling 
strength and initial state entanglement. For extended Werner states, we show that the concurrence 
is related to the difference of two functions: one is related to dephasing and the other longitudinal 
relaxation. A physical intepretation based on the generalized Bloch vector is given: revival only 
occurs for strongly-coupled noise and comes from the angular motion of the vector. 
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I. INTRODUCTION 

Entanglement is a property that sets quantum sys- 
tems apart from their classical counterparts [l[ . In recent 
years, it has drawn great attention as an important re- 
source for quantum information processing and commu- 
nication, such as quantum cryptography [2[, dense cod- 
ing Q , teleportation [1] and exponential speed-up of cer- 
tain computational tasks [U . Interactions with the noisy 
environment inevitably degrade quantum coherence and 
thus entanglement. It has been shown that although lo- 
cal (one-body) coherence decays continuously, global co- 
herence (entanglement) may terminate abruptly in a fi- 
nite time, a phenomenon known as entanglement sudden 
death (ESD) Q. To date ESD has been demonstrated in 
two different experiments Q. 

In the theoretical investigations of ESD, the Markov 
approximation is commonly used; namely, the environ- 
mental noise has short, or rather instantaneous, self- 
correlations pi-tll) ■ Non-Markovian noise, however, is 
widely observed in solid-state systems [HI, {l3j] and may 
even serve as the dominant source of decoherence jb4l - tl8l ] . 
It is thus crucial to extend the current understanding of 
ESD to the presence of non-Markovian environments. It 
is known that for noninteracting qubits coupled to their 
own independent environments, nonmonotonic time de- 
pendence of the entanglement can only occur with non- 
Markovian noise [19l - t2l| . Mazzola et.al. pointed out 
that a common structured non-Markovian reservoir pro- 
tracts the disentanglement process and enriches the re- 
vivial [22j . Yu et.al. considered a pure dephasing classi- 
cal Ornstein-Uhlenbeck noise model and found the short- 
time behavior of the entanglement evolution is markedly 
modified [23j] . Some progress has been made in modeling 
the behavior of interacting 2-qubit systems in the pres- 
ence of charge noise (24[. For the most part, however, 
the evolution of entanglement in a non-Markovian envi- 
ronment is still an open question. 

Pure dephasing noise models d, [2l], [23I [25I- 281 and 
variations of the Jaynes-Cummings model [1, 19, 20, [IH 
have been the standard testbeds for ESD. Our work 



however, is motivated by the recent discoveries in su- 
perconducting qubit designs flit - tltl ] where the working 
point of the qubit is tunable: an arbitrary mixture of 
dephasing and relaxational noise is possible. In these ar- 
chitectures, it is found that 1/f noise is the dominant 
source of decoherence. Thus a thorough understanding 
of those experiments requires models that a) deal with 
non-Markovian noise, especially random telegraph noise 
(RTN) [I3l which are the basic building blocks of 1/ f 
noise; b) treat dephasing together with relaxation [29j- 

n. 

In this paper we use the quasi-Hamiltonian method 
t29t - t3l| to investigate bipartite disentanglement of two 
independent qubits caused by uncorrelated sources of 
classical RTN. This method is extremely powerful for 
these types of problems and we will be able to obtain 
many analytic results, in an area of research dominated 
by numerical studies. Four issues are addressed. Firstly, 
it is known that RTNs can be put into two categories 
according to the ratio of their switching rates and cou- 
pling strengths to the qubit, namely the weakly-coupled 
(fast, Markovian) ones and strongly-coupled (slow, non- 
Markovian) ones [lo^tl^l • We thus seek for qualitative 
differences in the disentanglement caused by these two 
types of RTNs. Secondly, we exploit the working point of 
the qubit to see if this extra degree of freedom affects ESD 
and revival. This is particularly valuable when the work- 
ing point can be varied, since the number of strongly- 
and weakly-coupled noise sources can actually be tuned. 
Thirdly, some entangled states are known to be more ro- 
bust against disentaglement than others [25[. We thus 
examine different initial states, both pure (generalized 
Bell states) and mixed (extended Werner states). Fi- 
nally, we compare two noise models, the two-one model 
where only one qubit is subject to a RTN source and the 
two-two model where both qubits are attached to RTNs 
individually. This allows us to see the effect of noise lo- 
cality on entanglement, a global property. 

There are two distinct physical effects that lead to dis- 
entanglement. One is the movement of entangled states 
toward product states: this can happen even in the ab- 
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sence of noise. Second is the movement of entangled 
states towards completely mixed states: this requires 
noise. We will separate these two routes to disentangle- 
ment as far as we can by considering pure entangled ini- 
tial states and mixed entangled initial states; in the latter 
case, the distance from purity can be parametrized if the 
mixed states are chosen as generalized Werner states. 

We also propose the magnitude of the generalized 
Bloch vector |n| as an appropriate purity measure for 
multi-qubit states. In the single qubit case, states on 
the surface of the Bloch sphere (|n| = 1) are pure, while 
states at the origin (\n\ =0) are completely mixed. In the 
2-qubit case, the situation is not quite so simple, since the 
set of admissible states is not spherical. However, as we 
shall show, the pure states lie on the surface \n\ = \/3. 
In the 2-qubit state, we can then separate the sources 
of entanglement by computing both the concurrence and 
|n|. 

The paper is organized as follows. In Sec. HH we 
introduce the model Hamiltonian and apply the quasi- 
Hamiltonian method to reduce the two-qubit problem to 
the singlc-qubit problem. In Sec. IIIII we discuss the de- 
phasing and relaxation behavior of single qubit as func- 
tion of the qubit working point and RTN property. In 
Sec. IIVI we define the two-qubit entanglement measure. 
In Sec. [V] and I VII we solve the two-one and two-two 
models. In Sec. IVII1 we discuss the geometrical inter- 
pretation of the Bloch vector and qualitative difference 
between weak and strong coupling, and give overall con- 
clusions. 



II. MODEL 

The Hamiltonian of the system is given by 

H ^ = -\ E [Bov? + s K (t)g K -a K ]. (1) 

K=A,B 

where A, B refers to the two qubits, £?o is the energy 
splitting of the qubits between the ground and excited 
states, g is the coupling of the RTN to the qubit, and s(t) 
is the RTN sequence that switches between the values ±1 
with an average switching rate 7. Here a = [<j\ \ a%\ 0*3] is 
the triad of the Pauli matrices. 

This Hamiltonian is general enough to describe any 
qubit subject to c lassical R TNs. For a superconducting 
flux qubit, B = \J e 1 + A 2 , where e is proportional to the 
applied flux through the superconducting loop and A is 
the Josephson coupling, which are the energy difference 
and tunneling splitting between the two physical states. 
In this case 9 = tan _1 (A/e) is independent of K and the 
angle 9 is called the working point of the qubit. Since 
flux noise is along the e direction, 9 is the angle between 
the noise direction and the qubit eigenstate direction, 
and it can be varied by changing the applied flux. For 
spin qubits, 9k is simply the angle between the applied 
field and the magnetic noise field of the K-th fluctuator. 



For general problems, an arbitrary power and angular (9) 
spectrum can be obtained by superposing noise sources. 

In this paper we treat unbiased noise, so s(t) = 0, 
where the overbar denotes averaging over the noise dis- 
tribution. We also have 



s(t)s(t') =exp (-27 \t-t'\). (2) 

The noise autocorrelation function dies off exponen- 
tially, corresponding to a Lorcntzian power spectrum 
<Srtn(w) = 47<7 2 /(w 2 + 47 2 ). There are thus two time 
scales set by the RTN's characteristics: the correlation 
time of the environment r e ~ 1/7 and the noise induced 
looping time of the qubit 77 ~ 1/ (gcos9), where 9 is the 
angle between the energy axis z and the noise coupling 
direction g [3l[. The relative lengths of these two time 
scales are critical for the qubit decoherence and disentan- 
glement. Indeed, r e < ti is effectively the weak-coupling 
(Markovian) region while r e > Tg is the strong-coupling 
(non-Markovian) region. 

The density matrix pab (t) is 4 x 4 and can be expanded 
in the generators fa of SC/(4) 

Pab (t) = \(^i + Yl n * W ^ ' ( 3 ) 

where I4 is the 4x4 unit matrix. n,(i) = Tr [pAB{t)fa] 
are the components of the generalized Bloch vector ft; 
they are all real. The choice of the set of 15 Hcrmitian 
matrices {fa} is a choice of basis. We will take them to 
be 

Wa®<Jb} \ Wo® cr }, a, b = {0,1,2,3} (4) 

where 00 is the 2x2 identity matrix. We adopt the base-4 
ordering convention such that p\ = ao<g>ai, P2 = 0q®<J2, 
etc. The fa are an orthonormal basis for the density 
matrix space with the inner product (/x, /i') =Tr[/i/x'] /4. 
Note also that Tr fa =0. 

The rii (t) are not the most common way to character- 
ize a quantum state. However, they have a direct physical 
meaning. For example, since fao = (J2® &2, we have that 

"10 (t) = Tr \pwPAB(t)] = (02 ® er 2 ) . 

Thus nio is the value of a certain spin-spin correlation 
function. 

Furthermore, the rii (t) collectively form a measure of 
the purity of the state [34| . A pure state satisfies p 2 = p. 
Therefore any pure state satisfies =Tr (p — p 2 ) = | — 

j |n| 2 ,so \n\ = V3. At the other limit, the completely 
mixed state p = I4/A gives ft = and Tr(p — p 2 ) = 3/4. 
However, not all states with |n| < y/3 respect positivity 
34|. 

Using the quasi-Hamiltonian method [2^, [3(| , the time 
evolution of the quantum system in the presence of clas- 
sical noise can be cast into a time-dependent transfer 
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matrix T(t) acting on the generalized Bloch vector, 

n(t) = T{t) n(0) 

Note T(t) is real but is not orthogonal once the average 
over noise histories has been performed. Thus nonorthog- 
onality is a direct consequence of the incoherent environ- 
ment. 

To simplify notation, we include fj-o = I4 in the gener- 
alized Bloch vector, i.e. we define the 'extended' gener- 
alized Bloch vector as n = [no;n\. Note no(t) = 1 for all 
time. We thus have 

n(t) = T(t) n(0) 

and the time evolution T(t) can be succinctly written as 

T(t)=R A (t)®R B (t), 

where R K (t), K = A,B are the 4x4 'extended' single 
qubit transfer matrices. They are derived from the single 
qubit transfer matrices R K by padding the matrix: the 
(00) entry is set to 1 and the (0£), (£0) entries for i = 
1, 2, 3 are set to 0. Now 



R K (t) = (x f \exp(-iH«t) \i 



K = A,B 



and 



H. 



K 



-£7 + £7X1 



BqL z + T 3 gK ■ L 



Here and |x/) are related to the initial distributions 
of the RTN. We only consider unbiased RTN in this pa- 
per: \i f ) = \x f ) = [1; l]/y/2. 

If the qubit is not subject to any noise, the 'extended' 
single qubit transfer matrix can be written as 



1 0' 

cos B t sin B t 

— sin _B i cos B t 

1 



R (t) is orthogonal and thus conserves the length of the 
Bloch vector. The time evolution of the single qubit 
Bloch vector is simply a precession along the z direction 
with Larmor frequency Bq. 

If the qubit is subject to a RTN, the transfer matrix 
can be calculated by known methods [30|,[3l| and we have 



R 



RTN 



(*) = 



10 

( (t) cos B t C(t) sin B Q t 

-C(t) sin B t C(t) cos B t 

e 



(5) 



where ((t) characterizes the dephasing behavior of the 
qubit and Ti is the longitudinal relaxation rate. They 
will be discussed in the next section. i? RTN is generally 
a nonorthogonal matrix since both £ (t) and cxp(— Tit) 
decrease with time. 



As been pointed out by Bellomo et. al. [19| , the dynam- 
ics of -/V-qubit density matrix elements follows from the 
dynamics of each single qubit, and it is essentially inde- 
pendent of the initial condition of the total system. This 
feature is especially clear in our formalism, since T(t) is 
constructed directly from tensor products of single qubit 
transfer matrices R it) and is independent of the ini- 
tial conditions. It is know that collective channels (com- 
mon reservoir) can also lead to entanglement instability, 
either in the Markovian or non-Markovian environ- 
ment [22j. In this case, the two-qubit transfer matrix 
T_(t) is no longer of product form of single-qubit transfer 
matrices. The noise correlations glue up the single qubit 
Hilbert spaces and indirectly couple the two qubits [331 ■ 

The results in this paper have been calculated using 
the quasi-Hamiltonian method and have also been veri- 
fied through numerical simulations. A single numerical 
" run" is calculated by generating a sequence of RTN and 
then exactly solving the density matrix for that given se- 
quence. The final numerical simulation result is found by 
producing thousands of runs (40,000 for the figures in this 
paper) each with a different sequence of RTN, and then 
finding the average density matrix over all the runs. This 
allows us to numerical simulate the quasi-Hamiltonian 
results which are inherently averaged over all RTN se- 
quences. The numerical and quasi-Hamiltonian simula- 
tions are in agreement to within round-off error and could 
be plotted on the same graph without any visible differ- 
ence. 



III. SINGLE-QUBIT DEPHASING AND 
RELAXATION 



It is known that the fast (weakly-coupled, Markovian) 
RTNs and slow (strongly-coupled, non-Markovian) RTNs 
have qualitatively different effects on the single qubit 
time evolutions |30l - l33| . 

The two types of RTNs are separated by the criterion 
gcosO = 7 . For the fast ones (7 > gcosO), we have 77 > 
r e and the memory of the environment is short comparing 
to the looping time. Redfield theory (3(| applies in this 
case and both dephasing and relaxation of the elements 
of the density matrix is exponential at long times while 
the very short time behavior is quadratic [30(. 

For the non-Markovian RTN sources (7 < gcosO). 
T\ < T e , the correlation of noise is long enough to make 
possible the looping of the Bloch vector on the Bloch 
sphere. This looping manifests itself as oscillations in the 
Free-Induction signal (FID) and Spin-Echo (SE) signals 
P. 

The function £ (t) introduced in Eq(5]is directly related 
to the FID signal, i.e. npiD(i) = cosBot £(£). Physically, 
C (t) is the probability for the single-qubit Bloch vector 
to return to its starting point on the Bloch sphere in the 
rotating frame when no pulses are applied. We thus 
call ( (t) the dephasing function. Note ((t = 0) = 1 and 
^(t 00) — if dephasing occurs. 
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For the Markovian RTN we have the well-known re- 
sults: 



-r,t 



where 



C(t) = e 
T 1 g 2 cos 2 i 



r 2 = 



2 2 7 
2-/g 2 sin 2 



(6) 

(7) 

(8) 



4 7 2 + B 2 " 

For the non-Markovian RTN 

C(t) = e~ 7 * [cos(g cos 0t) + ei sin(5 cos 9t)\ (9) 

and the longitudinal relaxation rate is 

Ti = 2 7 e 2 sin 2 9, (10) 

where t\ = 7/3 cos 9 and e 2 = g/Bo. In each case the 
dephasing function f (0) = 1 and £ (t — > 00) = 0. 

Unlike the monotonic decay in the Markovian RTN 
case, £ (i) is oscillatory in the presence of non-Markovian 
RTN. It hits zero at discrete points in time 



C (*/) = at U 



7r^-tan _1 (l/ei) 
g cos 9 



1,2, 



(11) 



Note in both strong and weak coupling region, the lon- 
gitudinal qubit relaxation can always be well character- 
ized by a single exponential coefficient T\. 

Finally, we note that there is an exact result for £ (t) 
at the pure dephasing point 9 = [ITj] 



-Jt 



cosh ( v^7 2 — g 2 tj + 



sinh I \/j 2 — g 2 t 



l-( 7 ) 2 



(12) 



For EqfT2l if g > 7 (strong coupling), the hyper- 
bolic functions need to be replaced by the corresponding 
trigonometric functions 



C(*) = 



-fi 



cos 



(Vs 2- 7^) + 



sin 



I 1 1 



(13) 



Using the exact result, the zeros of £ (i) are given by 
tt£ — tan~ 



-, *=1,2,3,... (14) 



Vd 2 - 7 2 

We see Eq[TT]is indeed the the correct behavior as 7 — >• 0. 



IV. CONCURRENCE 

For bipartite entanglement, all entanglement measures 
are compatible and we use concurrence [13] for its ease 
of calculation. The concurrence varies from for the 
disentangled state to 1 for the maximally entangled state. 
It is defined as C AB (t) = max{0, q(t)}, and 



q(t) = Ai - A 2 - A 3 - A 4 



(15) 



where A; are the square roots of the eigenvalues of the 
matrix pabPab arranged in decreasing order and 

PAB = {<Jy 8> Oy)pA B {°t ^ CT * )> ( 16 ) 

where p* AB is the complex conjugate of the density matrix 

PAB(t). 

The product pabPab can be expanded as 



15 



PabPab — / 

i,3=0 



16 



where = (a y ® a y )p*(a y <£> <r y ). Note pi 



-Pi, if 



i = 1, 2, 3, 4, 8, 12 and = pi for other Vs. 

Thus our formalism allows us to investigate the time 
evolution of bipartite entanglement with the knowledge 
of the Bloch vector n(t). 

ESD occurs when q(t) < since the 'max' operation 
forces C AB = and C AB is not an analytic function. 
For the situations considered in this paper, we find that 
Ai and A2 have the same long time limit and A3 = A4. 
Thus ESD happens if lim^oo A3. i(t) =/= 0. In this work, 
we also find q(t) has a structure of \(\ — £, where £ is 
generally related to the initial state as well as the longi- 
tudinal relaxation process. Thus ESD happens whenever 
£(i — > 00) 7^ since always decays to zero in the 
presence of dephasing. 

We also compute (t) , the magnitude of the gener- 
alized Bloch vector. This is a measure of purity in the 
singlc-qubit case. In the two-qubit case it appears to 
track C AB to a large extent, and might serve as a poten- 
tial measure of both entanglement and purity for multi- 
ple qubit (> 2) systems when concurrence is no longer 
defined. 



V. TWO-ONE MODEL 

In this section, we consider the case where only one of 
the two qubits is subject to RTN. The 'extended' transfer- 
matrix is thus given by 



(17) 



In this model qubit B enjoys coherent time evolution 
and is stationary in the rotating frame. Thus all dis- 
entanglement and decoherence of the two-qubit system 
come from qubit A's interaction with RTN. We can think 



of the qubit B as forming a kind of reference frame for 
entanglement with qubit A. 
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Pure States 



We consider the generalized Bell states as initial states 



I*) 
I*) 



--a |00> + y3 111) 
=<*|01) + /3|10) 



(18) 
(19) 



where a is real positive, (3 = \\(3\\e lS and a 2 + ||/3|| 2 = 1. 
Note a = l/\/2, p = ±1/^2 gives the Bell bases 
and I**). 

With |$) or \^) as initial state, a complete analytic 
solution is possible. 

The time evolution of the generalized Bloch vector for 
initial state |$) is given by 



n 3 (t) =2a 2 - 1, 



(20) 



n 5 (i)=2af [cos(2B t)Rc /3 + sin(2B t)lm /3] , (21) 

n 6 (*) =2a C [cos(2B t)Im /? - sin(2_B *) Re 0\ , (22) 

no(t) =n 6 (i), (23) 

n w (t)=-n 5 {t), (24) 

ni 2 (t) =n 3 (t)e- rit , (25) 



n 15 (t) 



-r,t 



while the other rii (t) = 0. Here f (t) and Ti are given by 
EqU H]or[3 Uni depending on the coupling of the RTN. 

With 1^) as initial state, the non-zero components are 



n 3 (t) =l-2a 2 , 
n 5 (t) =2a Re /3 ((t), 
n 6 (t) =-2a Im/3 ( (t) . 
n 9 (t) = - n 6 (t), 

n w (t) =n 5 (t), 

n 12 (t) =-n 3 (t)e- Tl \ 

ni 5 (t) = 



-Txt 



(27) 
(28) 
(29) 
(30) 
(31) 
(32) 
(33) 



Note in this case the generalized Bloch is not dependent 
on Bo- This is because af + erf annihilates |^). 

For both initial states, the square roots of the eigen- 
values of pabPab are 



Ai = 
A 2 = 



i 



3 - rit + 2|C|) 
3 - rit -2|C|) 



A 3 =A 4 =« (1 _ e -r li) 



(34) 
(35) 
(36) 




Figure 1: Pure dephasing noise for the two-one model. Top 
panels (a) and (c): Square roots of eigenvalues of pp(t). Bot- 
tom panels (b) and (d): magnitude of Bloch vector n(t) above 
and concurrence C AB (t) below. These are calculated from 
Eqs. I20II37I and confirmed by numerical simulations. Both 
qubits are operated at the pure dephasing point and only 
one of them is connected to the RTN, i.e. gi = 0.1, 32 = 0, 
8 = 0. The initial state is set to In (a) and (b), 

7 = 0.005 and the RTN is in the strong-coupling region while 
in (c) and (d), 7 = 0.5 and the RTN is in the weak-coupling 
region. Ai is plotted as a solid blue line, A2 as a dashed green 
line; A3 = A4 = 0. n(t) is plotted as a solid red line and 
C AB (t) as a solid black line. Time is in the unit 1/Bo. These 
results are exact. 



(26) 

which gives the concurrence as 



C AB = 2oVl -a 2 max 0, | CC*> I 



1 - e-^ 



2a\/l-a 2 max{0, |C(*)| -£(*)} 



(37) 



The concurrence has a remarkable form. £ (t) > is 
only related to the longitudinal relaxation rate T% and 
varies between for the unrelaxed state and 1/2 for the 
fully relaxed state; it describes relaxation only, while C (t) 
is related only to dephasing. Entanglement exists only 
when |f I > £. The effects of dephasing (the decrease of 
|f |) and relaxation (the increase of f) are additive, and 
they race to disentangle the state. However, f (t) can be 
oscillatory, and then revival of entanglement is possible. 
Once the envelope of f (t) is less than f (t) , entanglement 
is gone for good. 

As a function of working point 9, Ti is hnite except for 
8 = 0: the pure dephasing point. At this point I\ = and 
f- (t) = 0. Two of the four eigenvalues vanish: A 3 ^ 4 = 0. 
If f (t) is a monotonic function (weak coupling), then 

C AB = 2oVl-a 2 |C(6> = 0)| > at all hnite times and 
ESD does not occur. This is seen in Fig. [T] (c,d). If 
f (t) is oscillatory (strong coupling) then ESD and re- 
vival occurs, as seen in Fig. fT] (a,b). In fact the revival 
happens an infinite number of times, since the envelope 
is exponential. We again note that the result at the pure 



dcphasing point is exact; see Eqs. [T2land[T3l 

For > 0, intermediate working points, Ti > and 
£(i — > oo) = 1/2. A3 = A4 increase monotonically with 
time. In fact all eigenvalues approach the same limit: 
^1,2,3,4 (t 00) = a\[\ — a 2 /2 as seen in Fig. [3](a,c). If 
6 > 0, then ESD is inevitable. All A's have a^/\-a 2 /2 
as their long time limit. At strong coupling, there are 
a finite number of revivals, while at weak coupling, ESD 
occurs at a finite time and no revival occurs. This generic 
behavior is the same at all intermediate working points. 

As we begin in a pure state, we have \n(t = 0)\ = 
a/3. At long times the magnitude of the Bloch vector 
approaches a finite limit in the dephasing case because 
n,3 = ((Tq (8> 03) , proportional to the expectation value of 
the z-component of the spin for qubit B, is independent 
of time. This is due to the cylindrical symmetry of the 
Hamiltonian: [H, erg <& 03] = 0. In the case 6 > 0, \n\ 
decays to zero, though this happens on a time scale longer 
than is shown in Fig. Q] Nevertheless, in both cases the 
time dependence of tracks the time dependence of 

C AB to a large extent. In particular, the oscillations in 
C AB (t) observed at strong coupling are also present in 
|n(£)|- These oscillations might naively be supposed to 
come from oscillations between entangled and product 
states; the fact that the oscillations are also present in 
\n\ means that the disentanglement is coming essentially 
from mixing, even though it is non-monotonic. 

Notice that \n\ , unlike C AB , is a continuous function 
of the elements of the density matrix p; as a result it does 
not suffer sudden death but rather decays exponentially 
The oscillations in C AB are not as long-lived as those in 
|n|. This is due to the fact that once the envelope of £ is 
less than £, ESD kills off C' AB . No such effect occurs for 
|n|. 



B. Mixed States 

We have found that for pure initial states, the disentan- 
glement comes mainly from mixing. Additional insight 
is gained by taking mixed initial states. We will use the 
"extended Werner states" 0, HH as initial states, 

wf(0) = r 1$) ($| + i_I/ 4 , (38) 

wf(0) = r |*><*| + i^Il 4) (39) 

where < r < 1 measures the purity of the initial states. 
r = 1 gives the (pure) generalized Bell states of the pre- 
vious section and r = is the fully mixed state p = I 4 /A. 
\n{t = 0)| = \/3r. 

The generalized Bloch vectors for the Werner states 
(n w &, n w *) are related to the those for the generalized 
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Figure 2: Mixed noise for the two-one model. Top panels 
(a) and (c): square roots of eigenvalues of pp(t). Bottom 
panels (b) and (d): magnitude of Bloch vector n(t) above and 
concurrence C AB (t) below. The curves are calculated from 
Eqs. I20II37I and confirmed by numerical simulations. Both 
qubits are operated at an intermediate working point and 
only one of them is connected to the RTN, i.e. g\ = 0.1, 
<?2 = 0, 9 — 7r/3, (j> — tt/2. The initial state is set to | ( 1 >+ ). In 
(a) and (b), 7 = 0.005 and the RTN is in the strong-coupling 
region while in (c) and (d), 7 = 0.5 and the RTN is in the 
weak-coupling region. Ai is plotted as a solid blue line, A2 as 
a dashed green line; A3 as a dottd red line, and A4 as a dotted 
black line. n(t) is plotted as a solid red line and C AB (t) as a 
solid black line. Time is in the unit 1/Bq. 



Bell states ( 7%\ tj>> , Tl\ ) by simple scaling 

ft w *(t) = r ff| $> (i), (40) 
n w f(t) = r n\ 9> (t). (41) 

For both initial states wf(0) and u>*(0), the A's are 
given by 

Ai=ra||0||(e+|C|), (42) 

A a =ra||/8|| (£- |C|) , (43) 
A 3 =A 4 = ra\\l3\\ £ (44) 

where 

^(1 + re- r i*) 2 - (1 + e- r i*) 2 r 2 (2a 2 - l) 2 
4rav 1 — a 2 

and the relaxation function is given by 

^(1 - re- r i*) 2 - (1 - e- r i*) 2 r 2 (2a 2 - l) 2 
Arayl — a 2 

which depends both on the relaxation rate ri and the 
initial conditions of the qubit. 
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Figure 3: Top panels (a) and (c) square roots of eigenvalues 
of pp(t). Bottom panels (b) and (d): Bloch vector n(t) above 
and concurrence C AB (t) below. The curves are calculated 
from Eqs. I40H45I and confirmed by numerical simulations. 
Both qubits are operated at pure dephasing point and only 
one of them is connected to the RTN, i.e. g\ = 0.1, g-z = 0, 

6 = 0. The initial state is set to p* with r = 0.5. In (a) and 
(b), 7 = 0.005 and the RTN is in the strong-coupling region 
while in (c) and (d), 7 = 0.5 and the RTN is in the weak- 
coupling region. Ai is plotted as solid blue line, A2 dashed 
green line, A3 dotted red line, A4 dotted black line, n(t) solid 
red line and C AB (t) solid black line. Note A3 = A4 7^ and 
ESD occurs in both cases. Time is in the unit 1/Bq. 



Figure 4: Top panels (a) and (c): square roots of eigenval- 
ues of pp(t). Bottom panels (b) and (d): magnitude of Bloch 
vector n(t) above and and concurrence C AB (t) below. The 
curves are calcuated from Eqs. 1401451 and confirmed by numer- 
ical simulations. Both qubits are operated at intermediate 
working point and only one of them is connected to the 
RTN, i.e. 31 = 0.1,g2 = 0, 9 = tt/3, <f> = tt/2. The initial 
state is set to pf + with r = 0.5. In (a) and (b), 7 = 0.005 
and the RTN is in the strong-coupling region while in (c) and 
(d), 7 = 0.5 and the RTN is in the weak-coupling region. Ai 
is plotted as solid blue line, A2 dashed green line, A3 dotted 
red line, A4 dotted black line, n(t) solid red line and C AB (t) 
solid black line. Note A3 = A4 ^ and ESD occurs in both 
cases. Time is in the unit 1/Bq. 



The concurrence is then 
C AB =i 



|0, 2roVl-a 2 (lC|-0} 



(45) 



If Ti ^ 0, i.e., if the qubit is operated at intermediate 
working point, the long time limit of £ is 



at 



00 



y/i -r 2 (2a 2 - l) 2 



4ra-\/l — a 2 



If Tx = (pure dephasing) 



(46) 



(47) 



Thus the only situation where ESD does not happen is 
when r = 1 and Ti = 0. ESD is thus essentially a uni- 
versal behavior if qubits are subjected to RTNs. 

In Fig. [3] and IH we can see in greater detail how the 
introduction of r, i.e. the interpolation with I4 in the 
initial density matrix, changes the situation. First, it lifts 
A3 and A4 to finite values. This causes ESD to happen 
even at the pure dephasing point, r is essentially a radial 
variable in n-space. Comparing Fig. [T] to Fig. [3] (pure 
dephasing noise), we find, in agreement with Eq. 1401 that 
ft is simply proportional to r. Thus the decay time of 
|n| is unchanged. The dependence of C AB on r, however, 
is more complicated. There is an overall proportionality, 



but the relaxation function £ also dpends on r. This 
means that the decay time of C AB is reduced when r = 
0.5 as compared to when r = 1 (pure initial state). This 
is in agreement with Eq. 25] Comparison of Figs. and 
U] (arbitrary working point) shows that this qualitative 
behavior does not change when we have mixed noise. 

This says something important about the geometry of 
the 15-dimensional space in which n lives. If we move 
radially in this space from a point with \n\ = v3 and 
C AB = 1 (a pure state with maximal entanglement) 
to the origin (the maximally mixed state) along a path 
determined by our model the entanglement diminishes 
monotonically, but not smoothly, to 0. If we move on the 
surface of the sphere of pure states (changing angular 
variables only), it is obviously possible to move contin- 
uously from an entangled pure state with C AB = 1 to 
a product (separable) state that would have C AB = 0. 
We may summarize this by saying that in our model pu- 
rity is a radial variable in n-spacc while entanglement 
is, roughly speaking, the product of the purity and an 
angular variable. 



VI. TWO- TWO MODEL 

In this section, both qubits are subject to RTNs and 
the two RTNs are not correlated and do not necessarily 
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have the same prameters g and 7. The transfer matrix 
for this model is 

T(t) = R A TN (t)®R B TN (t). 

We only consider the extended Werner states here since 
the generalized Bell states are included as special cases. 

With wf as initial state, the non-zero components of 
the generalized Bloch vector are 



n 3 (t) =r (2a 2 - 1) e 



-rf t 



(48) 



n 5 {t) =2r ( A ( B a [cos(2B t)Re + sm(2B t)Im 0\ 



(49) 



n 6 (i) =2r C C a [cos(2B t)lm - sm(2B t)Re P] 



n g (t) =n 6 (t), 
n w (t) = - n B (t), 

n 12 (t) =r {2a 2 - 1) e~ r ^, 

n 15 (t) =r e -( r "+ r ?>*. 



(50) 
(51) 
(52) 

(53) 
(54) 



With wf as initial state, the non-zero components of 
the generalized Bloch vector are 



n 3 (t) =r(l-2a 2 ) e^, 
n 5 (t) =2ra Re P ( A ( B , 

2ra Im ( A ( B , 
n 6 (t), 
nio(i) =n 5 (t), 

(t) =r(2a 2 - 1) e" 1 ^, 



n 6 {t) 
n g (t) 



ni2 
n 15 (t) 



(rf +rf )t 



(55) 
(56) 
(57) 
(58) 
(59) 

(60) 
(61) 



Note there is no B g dependence in this case. 

The square roots of eigenvalues of pabPab for both 
initial states are 



M=ra\\0\\ (i+\C A C B \) 

M = ra\\0\\ (C-ICVl) 
A 3 = A 4 = ra\\p\\ C 



(62) 

(63) 
(64) 



where 



(1 + re-V*+ T ? ) f ) 2 - r 2 (2a 2 - l) 2 (e" r ^' + e~ r f *) 2 



4ravl 



and the relaxation function is 



(65) 



I _ re -(rf+rf )ty _ r 2 (2a 2 _ 1)2 ( e -r*t _ e -rft 



B + \ 2 



4raVT __ a 2 



200 400 600 800 




Figure 5: Top panels (a) and (c): square roots of eigenval- 
ues of pp(t). Bottom panels (b) and (d): magnitude of Bloch 
vector n(t) above and and concurrence C AB (t) below. The 
curves are calcuated from Eqs. I48H67I and confirmed by nu- 
merical simulations. Both qubits are operated at pure de- 
phasing point and both of them are connected to the RTN, 
i.e. gi — 0.1, <?2 = 0.1, 9 = 0. The initial state is set 
to |$+). In (a) and (b), 7 = 0.005 and the RTN is in the 
strong-coupling region while in (c) and (d), 7 = 0.5 and the 
RTN is in the weak-coupling region. Ai is plotted as solid 
blue line, A2 dashed green line, A3 dotted red line, A4 dot- 
ted black line, n(t) solid red line and C AB (t) solid black line. 
Note A3 = A4 = and ESD does not occur in both cases. 
Time is in the unit l/£>o- 



The concurrence is given by 



AB 



C 



where 



max 



{0,2Wl-« 2 [|C AB (*)!-*(*)]} (67) 



C (t) = C A (t)C B (t). 

This simple product form is due to the independence be- 
tween the two qubits. 

The long time limit of the relaxation function is 



£(t 00) 



1 



Arayl 



if rf ^ o,rf ^ 0. 



(66) 



If T A = rf = 0, EqH7|i s recovered and r = 1 prevents 
ESD from happening, as seen in Fig. [5j 

The qualitative behavior of all quantities is rather sim- 
ilar in Figs. [UE]and[5l au referring to dephasing noise. 
This confirms a qualitative picture in which both qubits 
undergo a random walk in their respective Hilbert spaces; 
relative variables therefore also undergo a random walk, 
but faster. Quantities like C AB that depend on the rela- 
tive variables have a faster decay time. (Note the differ- 
ence in the time scales on Figs. Q]and[SJ) 

Finally we note the norm of the Bloch vector \n\ re- 
sembles concurrence in all cases considered. This can be 
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Figure 6: "Phase diagram" of the behavior of C AB (t), given 
'extended' Werner state as initial state. In the upper region, 
we have ESD and revival before C AB goes permanently to 
zero. In the lower region, C AB dies just once and for all. 



seen from the explicit expression 



M = 7 'a 



8a s (l-a 3 )«" fl ) 2 + e-''^ ^ ^ , + ^ ^•) 
\ +(1 



2a 2 



-2rf t 



-2rf t 



(68) 



Since the generalized Bloch vector fully describes the sys- 
tem, a geometric picture of n for entanglement might be 
possible. To our knowledge, such a description is not 
yet available except for some special parameterized states 
1391 or close to n = 1401 ■ 



VII. DISCUSSION AND CONCLUSION 

The most important conclusion of the paper is that 
the disentangling effect of non-Markovian, or strongly- 
coupled, noise and the effect of Markovian, or weakly- 
coupled, noise is qualitatively different. In Sec. [V] and 
Vll wc sec C AB {t) can take on two forms before ESD 
occurs: oscillatory and exponential, depending on the 
coupling of the RTN. By numerical exploration, we have 
constructed the "phase diagram" in Fig. [51 where the 
boundary is given by 



ah 



The oscillatory behavior arises from looping of n on 
the Bloch sphere. This can only occur if the noise is slow 
enough that the topology if the sphere is fully explored 
before the Bloch vector decays entirely. If the noise is 
fast, the relaxation moves n along a radial path to the 
origin no looping occurs. 

If the noise acts only on one qubit, the two-one model, 
the situation can be analyzed in some detail. Qubit B 
is not subject to RTN and is stationary in the rotating 
frame. It effectively serves as a reference and the two- 
qubit concurrence is fully determined by qubit A, as seen 
in Eq. E3 




Figure 7: Single-qubit Bloch representation of the two-one 
model. Initial states are the generalized Bell states. Inside 
the cone C AB = 0. 



In the single-qubit Bloch sphere picture, 



ICI 

,-rit 



Pa | sin0 A | 
p A cosO A , 



(69) 
(70) 



where pa and 9a are the length and polar angle of the 
three dimensional Bloch vector of qubit A. 

Given generalized Bell states as initial state, C AB = 
is equivalent to 



2/9a| sin( 



PA cos ( 



< 1. 



Geometrically, it means C AB = as long as qubit A's 
Bloch vector falls inside the cone shown in Fig. [7] 

In the two-two model, however, both qubits have non- 
trivial time evolution. This simple one-qubit picture for 
concurrence then does not work and one needs to treat 
the full 15-dimcnsional Bloch vector for the whole sys- 
tem. 

One important question is the relation of this work to 
previous results on single qubits (31rl33j . The oscillations 
that occur in \n\ and C AB are clearly related to the noise- 
induced looping on the single-qubit Bloch sphere. For 
example, they have the same period. However, in the 
single-qubit case these oscillations occur in the tails of an 
overall Gaussian decay. They are much more pronounced 
in |n| and C AB . 

Entanglement revival was introduced in Re f. IX II) and 
later on shown to exist in different systems |41|. It has 
sometimes been attributed to back-action from the non- 



Markovian environment 22 1. There is no back-action in 
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our model so this cannot be a general statement. Fur- 
thermore, the revival seen in the present work is a simple 
oscillation and does not arise from any constructive in- 
terference of multiple reservoirs. 

In conclusion, we used the quasi-Hamiltonian method 
to study the entanglement dynamics of two non- 
interacting qubits subject to uncorrelated RTNs, utiliz- 
ing the generalized (15-dimcnsional) Bloch vector n. This 
turns out to be very well suited to determining entangle- 
ment measures such as the concurrence, since n has a 
rather direct relation to entanglement. We found in our 
work that disentanglement caused by classical noise on 
2-qubit systems falls into two distinct categories. In the 
Markovian noise case familiar from perturbation (Red- 
field) theory, the motion of n is essentially radial and 
ESD happens except in special cases. The time scale of 
ESD is similar to the time scale of exponential decay of 
|n|. In the Markovian noise case there is a combination 
of radial and angular motion of n; \n\ typically shows os- 
cillatory behavior, while the concurrence undergoes ESD 



and revival. The quasi-Hamiltonian method provides us 
a flexible way to deal with independent qubits and un- 
correlated noises. The formalism for multiple qubits has 
been established in Ref. [30| and it is therefore straight- 
forward to extend the present work to this case. Other 
future work would be to explore the effects of inter-qubit 
coupling and noise correlations. 
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